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The s-wave scattering length for the Positronium-Positronium interaction is estimated semi- 
phenomenologically with the long-range van der Waals force and the short-range repulsive potential 
that represents the hard core between two positronium. The obtained value of the scattering length 
is a 0.44 nm, and its stability is also checked for different parametrizations. Using this value, the 
Gross-Pitaevskii equation can be fixed for the Positronium Bose-Einstein condensates (Ps BEG). The 
' static properties of Ps BEG are studied from the solutions of that equation. The phase-transition 

^\ \ temperature shift of the Ps BEG due to the Ps-Ps interaction is also evaluated with the mean-field 

^\ ■ approximation. 
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. I. INTRODUCTION 

Recent realizations of Bose-Einstein condensates for trapped Alkali-atom gas have made a great development in 
quantum physics of many-body system where the condensation phenomenon is one of the most important concepts and 
\^ . used in a vast area of physics such as atomic, laser, condensed-matter, nuclear, elementary-particle physics and also 
' in cosmology After the experimental success of the condensation, a lot of physical phenomena has been observed 
, and studied on Alkali-atom BEC that include collective excitations, sound propagations, quantum interference, etc 
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The success of the Alkali-atom BEC and the developments of the laser cooling and trapping techniques encouraged 
the study for the BEC of bose/fermi particles other than Alkali atoms: recently, the BEC of Rb isotope with the 
sympathetic cooling technique ||^ and the potassium fermi condensates have been observed [Q. The bose/fermi mixed 
condensates are also studied for the potassium isotopes [^,|| . 

The positronium (Ps), the bound state of the electron and the positron, has been also an interesting candidate for 
the BEC: 

pH ■ • The Ps is the lightest atom with the mass wpg ~ 2me where rUe is the electron mass: 10'^ times lighter than the 

Q \ hydrogen atom. It leads to the large value of the critical temperature Tq of the Bose-Einstein phase transition: 

O ■ in free-gas and continuous approximation, Tc becomes 10 times higher than that of hydrogen BEC for the 

trapping potential of same size {Tc (x m~^/^). The effect of the Ps-Ps interaction for Tc will be discussed in 
this paper. 

X 

' • The Ps has finite life-time and self-annihilates into photons: Ps 27 with Tppg ~ 0.125 ns for para-Ps 

^ and Ps 87 with T^pg ~ 142 ns for ortho-Ps. Thus, the Ps BEC should be an unstable condensate and 

decays with gamma-ray radiations, the coincident observation of which will give the precise measurements for 
position/momentum distributions of the BEC. Those radiated gamma-ray photons are also suggested to show 
new physical phenomena such as gamma-ray laser by spontaneous amplifications and the supperradiance. 

However, it should be noted that these very short life-times will make technical difficulties in experimental Ps BEC 
formations: rpg are in the comparable order with the relaxation time Trei ns in atomic phenomena so that the 
cooling methods that depend on the natural relaxation (like the evaporative cooling) might be less effective for the Ps 
BEC formation. At present, several faster cooling methods are planned by several experimental groups for Ps-BEC 
formations, and the preliminary experiments are progressing |8-10|. 



Originally, the Ps BEC was considered as e+-e~ condensates as a source of gamma-ray laser the e+-e~ 

plasma in the strong magnetic field at the neutron star surface and the storage ring |jl^ , and also as the source of 
the gamma-ray burst in astro-physics jw) . The Ps BEC itself has been suggested as a source of the gamma-ray laser 
[ p^ , but detailed studies for the properties of Ps BEC including the Ps-Ps interaction effects have not been done. 

In the binding structure, the Ps is similar with the hydrogen atom (H), where the main differences are mass of 
positively charged particle (mp/me+ ^ 2000) (and their stability). The BEC of poralized H has been reported with 
~ 10^ atoms |M. More similar object are excitons in semiconductors: the bound states of electron and hole in 
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valence band with mh/rrie ^ 0.7. The exciton has also finite life-time r 30 ns and decays into phonon, so that it 
is quite similar with the Ps. The BEC of excitons in CuO has also been reported [0. The experimental success of 
those BEC encourages strongly studies of Ps BEC. 

In this paper, we derive the Ps-Ps interaction potential with the long-range van der Waals interaction and the short- 
range cut-off that represents the hard core repulsion of positronium. The s-wave scattering length is estimated from 
the solution of the Schrodinger equation using this potential. With the given scattering length, the Gross-Pitaevski 
equation for the ground-state wave function (order parameter) of Ps BEC is given, and, solving this equation with the 
trapping potential, we study the ground-state properties of Ps BEC. We also estimate the Ps-Ps interaction effects for 
the critical temperature Tc in the mean-field approximation, which are shown to be negligible except large particle 
number BEC. 



II. EFFECTIVE PS-PS INTERACTION POTENTIAL 



i) van der Waals Potential for Ps-Ps Interaction 

We consider the system of two positronium Ps(l) and Ps(2). The spacial coordinates of the electron and positron 
composing Ps(l) are denoted by xi and Xq and those of Ps(2) are by X2 and x;,, with which the relative coordinates 
are defined by j = Xj — x^ . Especially, for relative coordinates between electron and positron, we denote 

ri = r^i, r2 = rb2- (1) 

The hamiltonian of the system is 

H = HcM + Hi + H2 + Hint, (2) 

where Hcm is the center-of-mass part that can be neglected in the derivation of the Ps-Ps interaction because of the 
translational invariance. The Hi_2 in (|^) are the internal parts of Ps(l,2) with the Coulomb potential: 

p? o ahc , , 

i?i,2 = — , (3) 

•m^ ri,2 

where ri 2 = |ri,2| and pi,2 are the momenta conjugate with ri^2- The a in (|3|) is a fine structure constant and 
a = e^/(47reo?ic) 1/137 in the SI unit system. The interaction part of Ps(l,2) is represented by 

H,nt =ahc(- + — -] , (4) 

^12 ra2 rbi J 

where — jr^ | and r = rab — |rab| is the distance between two positrons, Ps(l) and Ps(2). 

Let us consider the case where the distance between Ps(l,2) are very large: r = rab > ri,r2. In that case, the Ps-Ps 
interaction potential can be determined by the van der Waals potential ]l8| , p^ . 

Eq. can be expanded by 1/r, and we get the dipole-dipole interaction: 

H....--^{3^^^^-r.r.}. (5) 

When r ^ 0, the interaction term Hint ^ 1/r^ in (|^) can be evaluated with perturbation. The unperturbed part 
is just Ho = Hi + H2 defined by (^) and the zeroth-order unperturbed wave functions that diagonalize the Hq are 

V'nirn(l) 'S) V'n'/'m' (2) whcrc i/'nim(«) IS the cigenstatc of Ps(i), H.,%l)nim{i) = Eii^Tpnirnii), with the binding energy: 

E^^^ = -m,^o?l{\(n -f 1)2}. (6) 

The n, Z, m are principal, azimuthal and magnetic quantum numbers. In the Ps state function, we dropped the spin 
angular momentum part and spin quantum numbers S" = 0, 1 and S^, because the spin-spin/spin-orbital interactions 
are weak, and give a very small modifications that can be neglected in the present semi-qualitative estimations. 
From the second-order perturbation for Hint, we get the van der Waals potential: 

p(2) 2,. ^2 6 >^ |(0Ml|V-n;„(l))n(0|d2|V'»';'m'(2))P 

E^^ = -a (he) ^ , (7) 
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where di^2 = ^1,2 are dipole-moment operators of Ps(l,2) and |0) = l^/iooo)- Using eq. (||), the denominator of eq. ^ 
becomes 

4°^ + El^,^ - 4"' - 4'' - ?^[2 - (n + 1)-^ - in' + 1)-% (8) 

where uq = Tic/ (am^c^) ^ 0.0529 nm is the Bohr radius. The term [n + 1)^^ in the denominator of eq. (||) is very 
smah for large value of n, so that it can be neglected in the first approximation 

4°^+4?^-i^r-i^r-f^, (9) 

Using eq. (||) and the completeness condition: ^ ^ {^\di\ipnim) {4'nim\di\0) — (0|z^|0), eq. becomes 

£;(2) = -^^^^|(o|z^|o)p. (10) 

Using (0|z^|0) = (2ao)^ with the Coulomb wave function t/iqoo for the ground state, we get van der Waals potential: 

Ei^) = -6aM?^. (11) 



The error for the above approximation (|ll|) has been evaluated by numerical summation for for the H-H 
interaction When their results are rescaled for Ps potential, they obtained 6.47 instead of 6 in (1^). It suggests 
that the present results ( ITll ) should be a reasonable estimation. 

ii) Effective Ps-Ps interaction potential with the short-range hard core 

In the short range part, the real part of the Ps-Ps interaction should be characterized by repulsive hard core at 
the distance of atomic diameter where two positronium begin to overlap. In the present paper, we treat this hard 
core potential phenomenologically in two manners: A) the cut-off potential at r = rp, ^p^^pg (Fig- 2) and B) the 
Lenard-Jones potential with r~^^-hard core (Fig. 3): 



{2aof , ^(2ao)" 



-6^-^ + C 



(12) 



The phenomenological parameters rc and C in ^pg!^pg and Vp^ pg are fixed to reproduce the binding energy Bpg^ 
of the molecular state Ps2, which has not been confirmed experimentally but predicted theoretically /numerically with 
the hamiltonian (|^:^,^) . Recent results for Bpg^ by elaborate numerical calculations are around 0.435 eV (2l]-|23). 



To reproduce -Bpg^ = 0.435 eV with solving the Schrodinger equation for Vpg pg numerically, we obtain A) rc 
1.78ao (with the potential depth vq — ~6.13afic/ao) for ^pg^^pg, and B) C = 2.25 (uq = —2.0ahc/ao and r'c = 1.90ao: 
the distance where ^Ps_pg(''c) ^ 0) for ^g'^pg- It should be noticed that two phenomenological potentials have 
almost the same hard-core distances rc-,rQ ~ 1.8ao, which are consistent with the diameter of Ps, 2aQ. 

We should comment about the effect of Ps spin state. Because of the Pauli principle, the Ps2 bound state exists in 
singlet state, so that the effective Ps-Ps potential using the binding energy of Ps2 as an input parameter should be 
considered for singlet channel. However, because the rc of ^pg^^pg' is very close to 2ao that should be expected to be 
the hard core length for nonsinglet channel the present potential should be applicable also for these channels. 

III. S-WAVE SCATTERING LENGTH BY EFFECTIVE PS-PS INTERACTION 

For the Ps-Ps s-wave scattering, the relative wave function of two positronium becomes vE'(r) = rx{r) that satisfies 
the Schodinger equation: 

- ^^^M + ^Ps-PsWx(r) = Exir), (13) 

where /i is a reduced mass /i = TOpg/2 ~ rUe- 

The s-wave scattering length a is obtained from the asymptotic form of the E = solution for eq. ( p^ ) as 
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x(r oo) 1 . 

a 



(14) 



In fig. 4, we show the numerical results for the E — wave function x(^) of 6q. ( [l3|) with Vp^ defined in the 
previous section. The one peak at r ~ 3.5ao shows that (only) one weak bound state exists for this potential and 
it corresponds to the Ps2 molecular state whose binding energy has been used as an input for Vpg_pg. As seen in 
fig. 4, for r ^ 5ao, x(r) in fig. 4 are in the asymptotic linear region, so that the scattering length a can be read off by 
X{a) = 0: 

A) a = 0.440 nm for Vg^°pg, B) a = 0.437 nm for V"pg^_pg. (15) 

Both results are almost consistent, and it means that the s-wave scattering length (the low-energy limit) do not 
depend on the detailed form of the short-range potential. Existence of only one weakly bound state in Ps-Ps system 
is also essential for the rigidness of a given with the effective potential 

For the illustration of the effects by bound states for the scattering length, the s-wave scattering state function xi^) 
is shown in fig. 5 for the effective interaction potential fitted with H-H binding energy. That potential has the larger 
vo than the Ps-Ps potential and produces many excited molecular states that can be seen as the oscillating behavior 
of xl*") for r ^ bttQ. Its behavior depends on the short-range part of the potential given only in phenomenological 
manner in the effective potential, so that it leads to the impreciseness of the connected long-range behavior of x{i^) 
from which the scattering length a is given. For the Ps-Ps case, it should also be too rough approximation to treat 
the molecular bound state Ps2 itself as that of the effective potential given in this paper, but the long-range part of 
the scattering state should be more reliable because its short-range behavior is very simple and the potential has only 
one bound state that was used for the input parameter. The consistency between the scattering lengths given with 
^Ps^Ps ^Ps-Ps supports it. 

Currently-used s-wave scattering lengths for Alkali atoms are, 5.77 nm for ®^Rb |2^, 2.75 nm for ^^Na |^ and 
— 1.45 nm for '^Li [|6). The resuhant Ps-Ps s-wave scattering lengths in (|l5|) are a ~ 0.44 nm, which is smaller than 
those for Alkali atoms in one order but larger than that for polarized H atom (triplet), a — 0.065 nm It shows 
that the Ps-gas behaves more like perfect gas than the Alkali atom gas. The rough estimation for the Ps-Ps scattering 
length can also be made directly from the binding energy of Ps2 a ~ (2mei3pg^)~"'^/^ ~ 0.3 nm Q, which is also 
consistent with the present result. 



IV. PS BEC WITH PS-PS INTERACTION 



The dynamical behavior of BEC is represented by the ground-state wave function (order parameter) ■(/'(a;) and, for 
weak-interacting bose gas in the harmonic oscillator trapping potential at T = 0, it satisfies the Gross-Pitaevski (CP) 
equation ||^: 

'^it " " 2^^^'"^ + ^mps'^Hor^^ + S^IV'I'^, (16) 

where loho is the angular frequency of the harmonic oscillator potential, and N is the condensed bose particle number: 
the normalization condition for -0 is taken as 



d''x\iP{x)\^ = 1. (17) 

For low-density bose gas, the boson-boson interaction is dominated by the low-energy s-wave scattering, and the 
interaction constant g in (|l^) is given with the s-wave scattering length a Q: 

9= a. (18) 

mpg 

Using the dimensionless variables f — r/anoj i — ^Hot £ind — o^qV' with the harmonic oscillator length ano = 
^fij {m-pgUJHo), the GP equation (|l6|) becomes 



4 



where the dimensionless constant gN is 



gN = i^iV. (20) 
Thus, the solutions of GP equation are similar for the same value of gN. It should be noticed that the interaction 



effects come into eq. (19) through ^A^-term, so that they are more effective for large value of gN. 

In fig. 6, the solutions of eq. (hq) are shown with a = 0.437 nm in (|l|) and a; = 3.1 sec , for several boson numbers: 
N = 2.6x 10"^ (curve 2), 2.0 x 10^ (curv e 3), 7 .9 x 10* (curve 4), 2.5 x 10^ (curve 5). With a{Ps) = 44 nm for 
Ps, we obtain a{Ps)/aHO — 2.16 eV^^/^V^i^-HO, which is much smaller than a{Rb)/aHO — 7962 eV^^^^^/huJHO for 
Rb atom: this large difference comes not only from weak Ps-Ps interaction a{Ps)/a{Rh) ^ 1/130 but from the mass 
difference ^Jm-p^/mnb ~ 1/281. Because of the small a/ano for the Ps-Ps interaction, for large particle number 
A'' = 2.6 X 10'' (curve 2 in fig. 6), the calculated \il^{r)\ changes only ~ 16% from that for the ideal gas a — Q (curve 
1 in fig. 6). To realize 70% change for |V'(0)|, the Ps of = 2.5 x 10^ should condensate (curve 5 in fig. 6). For Rb 
BEG, because of large scattering length a = 5.77 nm, N = 6.7 x 10^ is enough for realization of the profile 5 in fig. 6. 

For the ideal gas with no interactions, the critical temperature Tq can be calculated in thermodynamical continuum 
limit §: 

kgT^ = Q.MhujHoN^'^, (21) 

where fcs is a Boltzmann constant. The particle interaction has an effect to shift the critical temperature: Tq = 
+ STc'. The STq can be evaluated with the mean-field and semiclassical approximation p9t : 

^ = -1.33^A^Va. (22) 

For Alkali atom BEG, eq. ( p2[ ) is consistent with experimentally observed temperature shift |^ (~ 2% shift for Rb 
BEG of iV - 10^). In the case of Ps BEG, eq. (|2|) becomes 

^ = -2.87eV-i/2y^iVi/^ (23) 

It gives very small shift for Tq for Ps BEG, and can be neglected in the first approximation except BEG of large N. 

It should be noted that the concept of Tc is exact for equilibrium (thermodynamical) system. For real experimental 
situations, the Ps BEG might be performed as nonequilibrium system due to the small Ps life-time that is almost 
the same order with relaxation time. In that case, the Tc evaluated here should be considered as estimation for the 
energy scale. 



V. SUMMARY AND DISCUSSIONS 



In the present paper, we studied the effective Ps-Ps interaction potentials with interpolating the long-range van der 
Waals interaction and the short-range repulsive hard core. The latter was taken phenomenologically in two kinds of 
parametrizations: A) sudden cut-off and B) Lenard-Jones (and cut-off) types. The Ps-Ps s-wave scattering length was 
calculated from the solutions of the Scrodinger equation with those effective potentials, and we obtained a ~ 0.44 nm. 
It is (2 ~ 10)-times smaller than that of the Alkali atoms, and larger than that of the polarized hydrogen atom. The 
stability of a for the different parametrizations of the short-range interaction was also checked. Using that value of the 
scattering length, we studied the ground-state wave functions for Ps BEG with solving the Gross-Pitaevski equation. 
The Ps BEG is shown to be more perfect-gas-like than Alkali-atom BEG because of weaker interaction. The critical 
temperature shift STc due to the Ps-Ps interaction was also estimated and shown to be very small because of the 
weak interaction and the small mass of Ps atom. 

As commented in the last section, the inclusion of the dissipative effects to the Gross-Pitaevski equation is essentially 
important for further studies of Ps BEG physics discussed in the introduction. It will be discussed in a future 
publication [|o|. 
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FIG. 1. Relative coordinates of e"*" and e~ in Ps-Ps system. 

FIG. 2. Effective Ps-Ps interaction potential Vp^p^ against r/ao- The ao = 0.053 nm is the Bohr radius. The rc is the 
cut-off length that represents the hard-core repulsive interaction. The vo in figure represents the potential depth. 

FIG. 3. Effective Ps-Ps interaction potential ^p^pg in unit of ahc/ao eV against r/ag. The ao = 0.053 nm is the Bohr radius. 
The short-range part is parametrized by the Lenard- Jones-type r~^^-potential. The cut-off r'(^ is defined by ^^Ps(^c) ~ 0, 
and the vo in figure represents the potential depth. 

FIG. 4. The profile of the s-wave Ps-Ps scattering solution x('") for l^p^pg in unit of ahc/ao eV against r/ao {ao = 0.053 nm). 
The s-wave scattering length a can be read off at xi^-) ~ 0, and a = 0.437 nm. For ^p^pg, we obtain almost the same result, 
and a = 0.440 nm. 

FIG. 5. The profile of the s-wave H-H scattering solution xi''') for ^p^pg against r/ao {ao = 0.053 nm). 
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FIG. 6. The space-distribution for Ps BEG at T = with the trapping harmonic osciUator potential {loho = 3.1 s ^) 

for several condensed Ps numbers: A'^ = 2.6 x lO'^ (curve 2), 2.0 x 10® (curve 3), 7.9 x 10** (curve 4), 2.5 X 10^ (curve 5). The 
curve 1 with gaussian profile corresponds to the BEG with no interaction. The normalization condition of ip is given by (p^. 
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